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Abstract 

One-dimensional Bose gases are considered, interacting either through the hard- 
core potentials or through the contact delta potentials. Interest in these gases gained 
momentum because of the recent experimental realization of quasi-one-dimensional 
Bose gases in traps with tightly confined radial motion, achieving the Tonks-Girardeau 
(TG) regime of strongly interacting atoms. For such gases the Fermi-Bose mapping 
of wavefunctions is applicable. The aim of the present communication is to give a 
brief survey of the problem and to demonstrate the generality of this mapping by 
emphasizing that: (i) It is valid for nonequilibrium wavefunctions, described by the 
time-dependent Schrodinger equation, not merely for stationary wavefunctions. (ii) It 
gives the whole spectrum of all excited states, not merely the ground state, (iii) It 
applies to the Lieb-Liniger gas with the contact interaction, not merely to the TG gas 
of impenetrable bosons. 
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1 Introduction 



Physics of ultracold Bose gases is a rapidly developing field of research due to recent remark- 
able achievements in experiment and intensive theoretical investigations (see reviews [1-4]). 
Among a number of important advancements in this field, one of the major steps forward 
has been the realization of quasi-one-dimensional Bose gases in elongated cylindrical traps 
and waveguides by tightly confining the transverse atomic motion [5-19]. The TG regime of 
strongly interacting bosons has been reached in a one- dimensional optical- lattice trap [20]. 
The one-dimensional TG gas of bosons moving freely was also recently reahzed [21] for ^'''Rb 
atoms by trapping them with a combination of two light traps. By changing the trap inten- 
sities it was possible to vary the ratio of the effective interaction to kinetic energy, achieving 
the TG regime, with this ratio reaching 5.5. The physical properties of low-dimensional Bose 
gases have been discussed in reviews [22,23]. However the important problem of the Fermi- 
Bose mapping in one-dimensional gases has not received the proper attention. It is the aim 
of the present short survey to compensate this deficiency by concentrating primarily on the 
Fermi-Bose mapping, first advanced over 40 years ago [24,25] for one- dimensional Bose gases. 
Particular emphasis will be placed on demonstrating the generality of this mapping, whose 
applicability is essentially wider than solely to the ground states of impenetrable bosons, 
as it is often considered in literature, when the mapping is reduced to an absolute-value 
relation. We shall stress that the mapping in its general form [24,25] is vahd for the total 
spectrum of excited states, for the time- dependent Schrodinger equation, and not only for 
impenetrable bosons, but also for the Lieb-Liniger gas with contact interactions of arbitrary 
strength. 



2 TG gas 

The mapping theorem was formulated for a quantum one-dimensional system of N bosons 
at zero temperature [24,25]. First, only the stationary states were considered. But, since 
the mapping procedure does not involve time, it was later stressed that the same mapping 
also applies for the time dependent many-body Schrodinger equation, which was employed 
for treating temporal interference properties of the one-dimensional hard-core Bose gas [26- 
31]. Therefore from the very beginning we may consider a time-dependent Hamiltonian 
H = H{xi, . . . , xjsf, t) where Xi e [0, L], with i = 1, 2, . . . , A^". The Hamiltonian is written as 
a general expression 

H^K + U+^Y.^iJ, (1) 
in which K is a kinetic-energy operator, having for nonrelativistic atoms of mass m the form 

the term U = U{xi, . . . , xj^f, t) describes any external potentials, such as confining potentials, 
generally allowing for the time dependence; and $jj = ^[xi — Xj) is a two-particle interaction 
potential. For bosonic atoms, the wave function 

^B = ^B{Xl-,...-,XN-,t) (3) 
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is symmetric with respect to the permutations of any Xi and Xj. The wave function (3) is a 
solution to the time-dependent Schrodinger equation 

i%^^ = H^. (4) 

In the stationary case, this reduces to the eigenvalue problem 

Hilj = Ei/j . (5) 

The two-particle interaction is assumed to contain a hard core of diameter a, which can be 
conveniently treated as a constraint on allowed wave functions 

ip{xi, . . . ,XN,t) = {\xi-Xj\<a) (6) 

for 1 < i < j < N, rather than as an infinite contribution to the interaction potential. 
Under constraint (6), the infinite hard-core potentials can be omitted from the Schrodinger 
equation, at the same time including all other possible finite interactions into the term U. 
Then Hamiltonian (1) can be contracted to the form 

H^K + U , (7) 

while imposing constraint (6) on the solution of the Schrodinger equation. The latter in a 
particular case of point-like impenetrable particles simplifies to 

^{Xi,...,XN,t)=0 {Xi = Xj), (8) 

where 1 < i < j < N. 

The problem of a one- dimensional hard-core gas was raised by Tonks [32], who considered 
the statistical mechanics of the classical high-temperature regime, while Girardeau [24,25] 
gave the solution for the quantum problem. That is why the one-dimensional system of 
impenetrable bosons is now commonly called the Tonks- Girardeau gas. 

The mapping theorem [24,25] can be formulated as follows. Let a wave function i/jp = 
ipF^Xi, . . . ,XN,t) be a solution to the time-dependent Schrodinger equation (4), with Hamil- 
tonian (7), possessing the fermionic antisymmetric property with respect to all permutations 
of any Xj and Xj, for i ^ j, and satisfying constraint (6). Introduce a unit antisymmetric 
function 

N 

A{xi, . . . ,xn) = '[[sgn{xi - Xj) , (9) 

i>j 

in which 

, . X { 1, a; > , 

'^"^^^ " R = 1 -1, x<o. 

Then the bosonic solution to Eq. (4) is given by the mapping 

. . . , xjv, = A{xi, ...,xn) iPf{xi, ...,XN,t). (10) 



3 



By this construction, function (10) satisfies the same hard-core constraint (6) as V'f- In 
the case of point particles, condition (8) holds automatically owing to the Pauli principle 
for ipp. pTinction (10) is totally symmetric under permutations of any Xj and Xj. It satisfies 
the initial and boundary conditions directly following from those for ipp. In the case of a 
stationary uniform system with periodic boundary conditions, the latter are preserved under 
mapping (10) if N is odd, but if is even, periodic (antiperiodic) boundary conditions on 
ipB require antiperiodic (periodic) boundary conditions on ipp. However for large ^ 1, the 
character of such boundary conditions becomes not important. The fermionic wave function 
ipF can be considered as corresponding to a fictitious system of spinless fermions, or better 
to say, to a system of real fermions with frozen spins aligned in the same direction. 

The Fermi-Bose mapping (10) is valid for the time-dependent wave functions, since the 
antisymmetric function (9) does not include time. The system Hamiltonian may contain any 
external fields and any other finite particle interactions in addition to the hard-core ones. 
For the stationary Schrodinger equation (5), the mapping apphes for the whole spectrum of 
all eigenstates. 



3 Ground state 



It is solely for the stationary ground state that mapping (10) reduces to a simplified form 

V'^(xi,...,Xjv) = |'0^(xi,...,Xjv)| . (11) 

As is evident, the absolute- value mapping (11) cannot apply to excited states or to the time- 
dependent case, when the wave functions are, generally, complex, whereas mapping (11) 
yields only real functions. In addition, for excited states (11) introduces unphysical cusps in 
ipQ arising from the requirement of orthogonality of different fermionic eigenstates i/jq , and 
positivity of (11) violates the requirement of orthogonality of different bosonic eigenstates 
ipQ. These defects are not present in the original mapping (10). 

Under the assumption that the only two-particle interaction is a zero-range hard-core 
repulsion, represented by the hard-core constraint (8), and there are no external potentials, 
the ground state can be found explicitly [24]. Since the fermionic wave functions vanish auto- 



matically whenever any Xi 



for i j, the constraint has no effect, and the corresponding 



fermionic ground state is that of the ideal gas of fermions, given by a Slater determinant of 
the lowest single-particle plane-wave orbitals. The exact bosonic ground state was found 
[24] to be a product 

1/2 



^0 {Xi,...,Xn) 



2N{N-1) 



in which ko plays the role of the Fermi wave vector, 

ko = np (p = 



N 

n 

i>j 



sm 



ko{xi — Xj) 



N 



N 

T 



;i2) 



(13) 



The ground-state energy can also be determined exactly [24,33,34]. For large N ^ 
can be easily obtained from the expression 

fh^e\ dk _ h^kl ^ 

ko V 2m / 27r 67rm ' 



En — L 



1, it 



(14) 
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with ko defined by the integral 

^^-^r-'i. (15) 



which yields ko from Eq. (13). Then the ground-state energy is 

£„ = ^iV. (16) 

bm 

The lowest excitations above the ground state have a phonon character [24] with the sound 
velocity 

c=^. (17) 
m 

The pair correlation function 

g{x,x') = L^I \i/j{x,x',X3,...,XN)\'^dx3...dxN, (18) 
Jo 

with the ground-state wave function (12), depends only on the difference x — x', so that 
g{x, x') = g{x — x'), and 

/ X sin^(A;oa;) , , 

For X < L, one finds [24] that 

/ N sin^(A;o2^) / x 

The vanishing of5'(0)=0atj; = reflects the hard-core nature of the two-particle interac- 
tions. 



4 Trapped gas 

The general mapping (10) holds true in the presence of any external potentials. The case 
of the harmonically trapped TG gas has been considered and an exact solution for the 
ground state has been obtained [31,35,36]. One-dimensional harmonic trap is described by 
the potential 

1 ^ 

(21) 

i=l 

The ground state of the Bose gas is given by mapping (11). The fermionic ground state 
is a Slater determinant of the lowest N single-particle eigenfunctions of the harmonic 
oscillator, 

_ exp(-a;V2Zg) /x 



where /q = yh/mw is the oscillator length and i?n(") is a Hermite polynomial. By rearranging 
the corresponding fermionic determinant, one gets [35] the Bose function 



(22) 
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in which the normahzation constant is 

2N{N-l)/2 1 

The single-particle density, normalized to iV, is 

p{x) = N j \'il){x,X2, ■ ■ ■,xn)\^ dx2 . ..dxN , (23) 
which for the ground state (22) gives 

AT-l 

pix)^^Mx)\'. (24) 

n=0 

The pair correlation function becomes 



Af-1 



where 

A{X,X') = ^ ^l{x)^Pn{x') . 
n=0 

For large iV 3> 1, one has A(a;, x') ~ 5{x — x'). But for x = x', as is evident, 

A{x,x) = p{x) , 

because of which function (25) vanishes, g{x, x) — 0, as it must be for impenetrable particles. 

5 Type of order 

What type of order exists in the system is characterized by the behaviour of reduced density 
matrices [37-39] . Of particular importance is the first-order density matrix 

pi{x,x') = ^ j i'{x,X2, . . . ,XN)i'*{x',X2, . . . ,Xn) dx2 ■ ■ -dxN , (26) 

normalized as 

I pi{x,x) dx = N . (27) 

The presence or absence of long-range order is described by the properties of the eigen- 
values rij of the density matrix (26), which are given by the equation 



pi{x,x')(pj{x') dx' = njipj{x) . (28) 

The corresponding eigenfunctions ^Pj{x) are called the natural orbitals [39], since in their 
terms the single-particle density matrix acquires a diagonal representation 

Pi{x,x') = ■ (29) 
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The eigenvalues rij play the role of the occupation numbers of the related orbitals and are 
normalized as 

3 

which results from normalization (27). The Fourier transform of matrix (26) gives the 
momentum distribution 



i{k) = J pi{x,x') e-''^^^-^') dxdx' , (30) 

with the normalization 



/ 



n{k) - = N 

ZTT 



The problem of calculating the first-order density matrix (26) for a uniform gas was first 
considered by Schultz [40], who found it in the form of a Toeplitz determinant. Using the 
known asymptotics of the Toeplitz determinants, it was possible to prove the absence of 
Bosc-Einstein condensate by showing the power-law decay of the density matrix at large 
distance, that is, by demonstrating the absence of long-range order. The precise form of 
this power-law decay was found later by Lenard [41,42], who obtained the long-distance 
behaviour as 

pi(a;,0)~C-^ (|a:|^oo), (31) 

where /cq = vrp. The coefficient C = 0.92418 was found by Vaidya and Tracy [43,44]. Higher- 
order terms in the asymptotic behaviour of pi{x, 0) were derived by Jimbo et al. [45]. The 
most accurate results are due to the recent work by Gangardt [46], who obtained 



pi(x,0)~ 



y/koX 



cos{2kox) 3sin(2A;oa;) 



32(A;ox)2 8(A;ox)2 16(A;ox)= 

33 93cos(2/,-o.x-) 



+ 



2048(A;o2;)^ ' 2bQ{kQxY J ' ^^'^^ 

He also found the finite-size corrections for atoms in the harmonic trapping potential and 
for the case of circular geometry [46] . 

The investigation of the properties of the first-order density matrix (26) revealed that the 
number of particles in the Bose condensate, which is associated with the largest eigenvalue 
of eigenproblem (28), is of the order of A^o = sup^ rij ~ V^. This imphes that there is no 
real Bose-Einstein condensate in the uniform Tonks-Girardeau gas. 

The case of trapped atoms does not allow for a simple analytical expression of the largest 
eigenvalue rij. The multidimensional integral (26) was evaluated numerically by Monte Car- 
los integration [35,36]. Highly accurate results for large values of were found by Forrester 
et al. [47]. These results show that again, as in the spatially uniform case, A^'o ~ V^. 
Thus, there is no true Bose-Einstein condensate in the trapped TG gas, because of which 
the Gross-Pitaevskii equation, presuming the existence of a well-defined order parameter as- 
sociated with genuine Bose-Einstein condensate, has limited utility, especially, for temporal 
processes [26,31,48]. This is contrary to the case of the trapped ideal gas, where Bose-Einstein 
condensation can develop, though as a gradual crossover but not as a sharp phase transition 
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[49]. Nevertheless, since for iV 3> 1 the momentum distribution of the Tonks- Girardeau gas 
exhibits a peak n{k) ~ k~^^'^ in the neighbourhood of zero momentum, and a kind of order 
does exist, such a system displays some coherence effects, and the Gross-Pitaevskii equation 
has a limited region of applicability [50] . 

In order to accurately classify the type of order arising in the TG gas, it is possible to 
resort to the notion of the order indices, introduced for density matrices [39] and generalized 
for the case of arbitrary operators [51]. The order index of an operator A is defined [51] as 



log \TrA\ 

where || ■ || means a Hermitian norm. For an n-th order boson density matrix pn, one has 
||Pn|| ~ IIpiII". Also, log|Trp„| ~ nlogN, when N ^ 1. Therefore, applying definition (33) 
for a reduced density matrix p„, we have 

'^(Pn) - ■ (34) 



Taking into account that 
Eq. (34) can be rewritten as 



logA^ 

IPiII = maxrij = A^o , 
j 



log AT 

For the Tonks- Girardeau gas, both uniform as well as trapped, A^o ~ . Thence 

^{Pn) = ^ • (35) 

In the case of a genuine Bose-Einstein condensate with long-range order, one would have 
^(Pn) = 1- The order index (35) characterizes a system with mid-range order [39,51]. The 
occurrence of mid-range order means that, though there is no true Bose-Einstein condensate, 
some partial coherence does exist in the system. 



6 Lieb-Liniger gas 



Up to now, the Fermi-Bose mapping (10) has been applied to the TG gas, that is, the one- 
dimensional gas of impenetrable bosons. However, as it turned out, the applicability of this 
mapping is much wider, being valid for a large class of one-dimensional systems called the 
Lieb-Liniger gas, which is characterized by the contact two-particle interaction 



$(x) = ^o6{x) 



(36) 



The one-dimensional system with this interaction was studied by Licb and Liniger [52,53]. 
The delta potential (36) leaves the wave function continuous but yields a jump in the deriva- 
tive according to the condition 



d_\ ^ 

dxj dx-i I 



=a;,+0 



_d_ 

dxi 



IT 



(37) 
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The dimensionless ground-state energy of this gas, 



expressed as a function of the dimensionless couphng parameter 



e{g) = -zTTT-r ' (38) 



9 = ^. (39) 

is given [52] by the equation 

<9) = [f) jj{x)x'dx, (40) 

in which the function f{x) satisfies the integral equation 

2./W = 1 + 2a£^;4(^, (41) 

and the constant A is defined by the normalization condition 

I £ f{x) dx^l. (42) 

Numerical solution for e{g) was given in Refs. [52,54]. A detailed table can be found on the 
website [55]. An analytical asymptotic expansion in the weak-coupling limit reads as 

e{g) ~ y + 03^=^/' + C4^' + C5^'/' . (43) 

The coefficients 

4 1 29 

C3 = = -0.424413 . C4 = — ^ = 0.065352 

37r 27r"' 

were found by Lee [56,57]. The coefficient C5 is not known exactly. Its estimate is C5 = 

—0.017201. The strong- coupling expansion, being based on the numerical results [55], can 

be derived as 

, ^(. 4 12 32 80\ 
e{g) ~ e(oo) - z + ^ - ^ + -A ^ (44) 



where 



g g" g^ 



e(oo) = - (45) 

is the TG limit [24]. The first-order density matrix and the pair correlation function were 
investigated by Monte Carlo techniques [58]. It has been mentioned that mapping (10) 
can serve as a reasonable approximation for limited time intervals of the time-dependent 
Schrodinger equation [59]. 

Cheon and Shigehara [60,61] showed that mapping (10) is exact for the Lieb-Liniger gas, 
provided that the fermionic wave function satisfies the condition 

=^^. =-^^. . (46) 
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Here the derivative is continuous but the wave function is discontinuous. 

Recently a new development has been proposed [35,62] for the case of the spin-aligned 
Fermi gas, suggesting exploitation of the generalized Fermi-Bose mapping [60-64] in the 
opposite direction, by mapping the fermionic Tonks- Girardeau gas, a spin-aligned Fermi gas 
with strong one-dimensional atomic interactions mediated by a three-dimensional p-wave 
Feshbach resonance, to the trapped ideal Bose gas. 

It can be shown [62,64] that in the low-energy domain the one-dimensional longitudinal 
scattering of two spin-aligned fermions, confined in a single-mode harmonic waveguide, can 
be well represented by the contact condition 



F\xi=Xj+0 



aiD 



dipt 



dx 



(47) 



— *C J —1—0 

Here the following notation is used for the effective one-dimensional scattering length: 



aiD 



1 + 



3C(3/2) fa^^^' 



2V27r V/. 



n -1 



(48) 



in which is the p-wave scattering length, C(3/2) = 2.612 is the Riemann zeta function, 
and Z_L = ^jTi/muj±_ is the transverse oscillator length. The scattering length am diverges at 
{ap/l±Y — —1-134. The fermionic TG gas regime occurs in the neighbourhood of this reso- 
nance. The one-dimensional scattering lengths are invariant under the Fermi-Bose mapping 
(10), as a result of which the scattering length am is the same for bosons and fermions. 
The contact condition (47) is generated by the one-dimensional pseudopotential operator 

[62] 



$^(a;) 



1 (dip! 



dx 



x+O 



dx 



(49) 



x-Q^ 



in which 5'{x) is the derivative of the Dirac delta- function and the effective coupling strength 
is 



2h 



2 <^1D 



m 



This should be compared with the bosonic interaction strength 

• 

maio 

The spin-aligned Fermi gas maps to the Lieb-Liniger Bose gas, with the fermionic and bosonic 
interaction strengths inversely related [60,61] as 



4^ 



(50) 



From here, it is clear that a strongly interacting Fermi gas can be mapped to a weakly 
interacting Bose gas. In the limiting case, when at the resonance the fermion interaction 
becomes divergent so that am — > — oo, the corresponding Bose gas is asymptotically free. 
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7 Discussion 



In the present survey, we have considered one- dimensional Bose gases. The properties of 
such gases are drastically different from those of their three-dimensional counterparts. For 
comparison, we may recall the ground-state energy of the three-dimensional dilute Bose gas 
with a hard-sphere interaction. The dimensionless ground-state energy is 

2mEo 



~ 47ra^/^ (l + bia^/^ + 62" + b'^a In aj , 

where a = pa^ <C 1, with a being the sphere diameter, and where the coefficients bi and 63 
are 

61 = = 4.814418 , 6' = 8 f— - v^") = 19.653915 . 

ISa/tt V 3 / 

The coefficient bi was found in Refs. [65-69] and 62 in Refs. [70,71]. The coefficient 62 
has not been determined exactly. According to Wu [70], one has 62 = 62ln(127r) = 71.337. 
Hugenholtz and Pines [72] give 62 ~ 74.617. Another quantity that could be compared with 
its one-dimensional counterpart is the pair correlation function of the three-dimensional 
dilute gas of hard-sphere bosons, which is found [67] to be 



where r ^ oq, a <S 1 and 



ttQ being the mean interparticle distance. The behaviour of both Eq and g{r) for the three- 
dimensional gas is essentially different from their onc-dimensional analogs. Being principally 
different, the three-dimensional gas does not allow for a direct mapping between bosons and 
fermions, as in the one-dimensional case, which makes its treatment much more involved. 
While in the one-dimensional case, the Fermi-Bose mapping (10) gives a great advantage in 
the simplification of the mathematical description. 

Nowadays one-dimensional Bose gases are not just artificial models, but vise versa are 
real physical objects realized in a number of experiments [5-21] with quasi-one-dimensional 
traps and waveguides. Quasi-one-dimensional dilute Bose gases can be treated as purely 
one-dimensional gases, with their effective interaction of the contact type (36), where the 
interaction strength 



2ra 



1 



IC(l/2)|as 
V2l± 



is expressed through the three-dimensional scattering length and the oscillator length 
l± = ^h/mu!± of the confining radial potential [73,74]. Here C(l/2) = —1.4603. 

An additional dimension is provided by the Feshbach resonance techniques which make 
it possible: to vary the scattering length in a very wide range, transforming the interaction 
strength from weak to strong coupling; to realize heteronuclcar resonances between two 
different atomic species [75,76]; and to perform a crossover between a degenerate fermionic 
gas and a gas of bosonic molecules (see Refs. [23,77-79]). In quasi-one-dimensional traps, the 
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Feshbach resonance techniques would allow for a continuous tuning of the system properties 
from a weakly interacting gas to the hard-core TG gas [80]. 

The remarkable generality of the Fcrmi-Bosc mapping (10) provides us with a convenient 
practical tool for considering different regimes of one-dimensional or quasi-one-dimensional 
dilute Bose systems, from the TG gas of impenetrable bosons to the Lieb-Liniger gas with 
contact interactions. Moreover, mapping (10) is applicable for describing nonequilibrium 
processes in such one-dimensional gases. A great variety of available experiments, to which 
the general mapping (10) is applicable, makes the latter of fundamental importance. 
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